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Abstract 

The investigation on the production of particles in slowly varying but 
extremely intense magnetic field in extended to the case of axions. The 
motivation is, as for some previously considered cases, the possibility that 
such kind of magnetic field may exist around very compact astrophysical 
objects. 
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1. Statement of the problem 

A magnetic field of huge strength can give rise to real particles even if 
its rate of variation is very small: this posibility could be of some interest 
from a pure theoretical point of view, but it gains more physical relevance 
if one accepts that such kind of field configurations may be present around 
some very compact astrophysical objects [1-3]. In this case the time vari- 
ation is related to the evolution of the source, by collapse, rotations or 
else, and it is therefore very slow, in comparison with the times typical 
of elementary-particle processes. We can call the former time the macro- 
scopic time and the latter the microscopic one. The production of light 
particles in these processes has been analyzed in some detail in some pre- 
vious papers [4], with the suggestion that it is one of the mechanisms at 
work in the phenomenon of gamma-ray bursts[5]; the typical microscopic 
time is related to the electron mass since photons are produced through 
real or virtual intermediate states of e~ — e + -pairs. The lightest particles 
that could be produced are massive neutrinos, but the magnetic-moment 
coupling induced by the standard electroweak interactions is extremely 
small. 

There is, at least in the theoretical realm, another very light parti- 
cle, that is the axion[6]: owing to its dynamical characteristics it must 
be coupled also to the electromagnetic field [7], even more its electomag- 
netic coupling is being actively studied from an experimental side [8] and 
the possibility of detecting such particles as coming from nonterrestrial 
sources has already been foreseen [9]. It is immediately seen that the pro- 
duction of axions by a varying magnetic field must be realized through a 
mechanism different from the previously considered one, in fact the axions 
are coupled [7] only to the pseudoscalar density E • B so the presence of 
an electric field is necessary as a starting point, but a nonstatic magnetic 
field creates always an electric field and even though the rate of variation 
is small, the very large magnetic strength makes the electric field not a 
tiny one. 

In the present paper the coupling of the axions field with a given 
E ■ B density is written in standard second-quantized formalism, then the 
effect of time variation of that density on the axion vacuum is determined 
and the consequent production is calculated. The result depends both 
on the spatial shape and on the time variation of the magnetic field: in 
accordance with the prevailing astrophysical hypotheses [1-3] the magnetic 
field is seen as a bundle of lines of force which may safely be considered 
straight in comparison with the microscopic scale. The time-variation 
could affect both the shape and the strength of the fields, both are effective 
in the production process. The calculation procedure is not the standard 
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adiabatic approximation [10] as used in previous investigations [4], but the 
feature that one has to deal with a two-scale problem is still fully relevant. 



2. General form of the production probability 



The starting point is a second-quantized axion field in presence of 
given, classical, magnetic and electric fields. The axion field 4>(x) is cou- 
pled to the pseudoscalar density * G(x) = E(x) ■ B(x) and the coupling 
constant, of dimension of length, is here indicated by C. We assume that 
the interaction lasts from an initial time t Q until a final time t. 

So we have for the axion field the expression: 

<P(x) = (j) {x) + X (x) =M%)+C J A R (x - y)G(y)d A y (1) 

Here Ar(x — y) is the standard retarded Green function, the source is a 
c— number, so the same holds for x- The field is free before t D , where it 
has the standard expansion: 

1 r rfik r i 
M*) = -—^ / 4= a (k)e ik - r -^ + at (k)e- ik ' r+i ^ , (2) 



(2tt)3/2 7 ^J~ k 

then it acquires a contribution from x, this term has the following actual 
expression: 

*(*> = (S^TS / g I 8 """" £e'"^M*e"" - ex.] (3.1) 

' Mt)= (2^572 /cMle-'"- (3.2) 

The reality condition for G, that gives gl = g-k have been used, together 
with the initial condition x(t ) = 0. The separation into positive and 
negative frequencies in x( x ) is unambiguous until the typical frequencies 
of G are small. 



Having found the time evolution of the field the total production of 
axions is calculated in Heisenberg description of motion, i.e. we take an 
initial state a Q (k)|o >= Vk, i.e. the vacuum in the absence of inter- 
action then we express the mean particle number as the time dependent 
term 

Af(M) =< o\a\\t,t) a(k,£)|o > 

* the same kind of coupling is possible also for the neutral pion, but in this case 
other channels of production are present [11] 
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Since all the effect of the interaction is a c— number shift on the opera- 
tors a(k, t) = a D (k) + 6(k) the calculation is easy, in particular when the 
interaction no longer acts we get 

AT / (k) = |6 / (k)| 2 (4) 

The expression of 6/(k) can be read off from eq.s (2,3.1,3.2), it is: 



^ k) = vfef el "^ (T) * 



(5) 



3. Detailed calculations 

As anticipated in the introduction a more definite model of the mag- 
netic field can by a field of uniform direction (at a given time) with some 
transverse shape, both the direction and the shape may vary in time, one 
possible restriction is the conservation of the total flux. More explicitly 
these conditions are realized by giving: 

F ii A v F 

A=- o—[l-w(iJ,r±)] B = n d±w(fxr±) . (6) 

z7T r]_ 7rrj_ 

The unit vector n gives the instantaneous direction of the magnetic field, 
r± = \lr 2 — n • r 2 and d±_ the corresponding derivative, F is the total 
magnetic flux; the parameter p defines the size of the field in the transverse 
directions, w is taken to be cylindrically symmetric and, obviously it must 
go to zero at infinity, the requirement w(0) = 1 avoids singularities in E. 

Since E = — A and A ■ B = only the terms coming from the time 
variation of the direction of A contributes to G = E • B so we get: 

G=Q 2 ^[(l-v,)d ± w] (7) 

Here J gives the angular velocity of n i.e. J = n A n. In this configuration 
the Fourier transform of the source is 

9k(t) = ^^C(2F) 2 5(n ■ k) J ■ kS(k ± /ri (8.1) 

S(k ± /n) = (jJL/k ± ) J J^pkJriKl-wip^w'ip^dp/p (8.2) 

Here J\ is the Bessel function of order one and w' indicates the derivative 
with respect to the argument. It is useful to remember that, owing to the 
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presence of the factor 8(n ■ k) in the expression of S we can substitute k\ 
simply with k 2 . 

We now remember that the model of magnetic field we have at hand 
is such that it is uniform along one direction, but this direction is contin- 
uously varying, so a most significant quantity is obtained by an angular 
integration 

|£ = k 2 J dn k Af f (k) (9) 

So we need a quantity like J dQkgk(T)gl(T') which contains a singularity 
due to the presence, in the domain of integration, of a 5-square term which 
arises for r = r'. This is clearly due to the unphysical assumption that 
at every time there is a direction in which the pseudoscalar density G 
is absolutely uniform. Since we are integrating over the direction at the 
end the effect of the singularity is mild enough, it results in a logarithmic 
divergence, however this fact must be explicitly death with, considering a 
finite extension of the fields. A more careful treatment is mathemetically 
heavier, so it is presented in an Appendix. 
The final result is given by 

Ft A k 2 

^ ^C 2 (2F) 4 ^-^[S(k/ij)] 2 V2ln4kL. (10) 



If we give a definite transverse shape to the fields we get, evidently, 
a definite answer. In the situation where the fields change direction but 
not intensities, i.e. we keep a constant, we may give, tentatively the form 

2 

w(p) = e~ p . Then the expression for the function S(k/fj) is[12] : 



S(k/fi) = 




The limit k — > of this expression is finite, so the whole production goes 
to zero only owing to the phase-space factor k 2 in front of the expression in 
eq.(10). This result, as it appears from the whole derivation, can be valid 
only for axion masses and so energies which are definitely larger than the 
typical frequencies of the astrophysical phenomena, no resonant dynamics 
is included. 

4. Some conclusions 

The rate of production of axions by a slowly varying but very strong 
magnetic field has been calculated. The conditions are such that the as- 
trophysical frequencies are very much lower than the proper frequencies of 
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the axion field. In fact with an axion mass tua of the order of leV[13] the 
typical frequencies of the field shall exceed 10 15 Hz and so no resonance 
conditions appear realistic. The relevant parameter turns out to be the 
dimension of the spatial inhomogeneity, in the present model here 
also is seems very reasonable to assume that /j « rriA- It is also possible 
to give an expression for the total number of produced particles. 

f 2 F 4 // 3 

A = C~ — VlnL/j (] 

m A 

Some of the factors owe their origin to the general form of the interaction, 
eq(l) and to dimensional requirament, it appers clear the role of the total 
rotation of the field (\&) in determining the overall production; so when 
the rotation is uniform the rate is proportional to the angular velocity. 

The numerical factor is more model dependent, in the chosen case it 
is C = 80^2^ - 2 - V2] = 0.704 . . .. 

The transformation from the incoming Heisenberg field to the outgo- 
ing field can be implemented by the simple unitary operator 



U = exp 



d 3 k[a(k)b*(k) -a f (k)6(k)]] 



as: 



Ua{k)U ] = a(k) + b(k) . 



The actual form of the evolution operator gives the further information 
that the axion are produced with a Poissonian distribution of multiplicity, 
strictly speaking this is true for a production in a totally defined state, 
in operations like the one leading to eq. (9) this particular form can be 
blurred. 



A very short mention on this problem was made at the XXXVI In- 
ternational Symposium on Multiparticle Dynamics, Paraty, R.J. - Brazil, 
September 2006 



Appendix 

We want to calculate / dQkgk(T)gl(T'), with the functions g given by 
eq.(8.1) and taking care of the finite extension of the magnetic field. This 
is implemented by substituting the 5— functions as: 

S(n-k) ^Lyr-^expt-LV-k) 2 ] (A 

The calculation is performed in the particular case in which the rotation 
takes place in a constant plane, so J|| J'. The integration over the angles is 



6 



(A.2) 



performed in Cartesian coordinates, it is useful to introduce the unit vector 
of the three-momentum direction k = kv, then J dQk = 25(t> 2 — l)<i 3 i> 
and through standard although lengthy calculations the representation is 
obtained: 

J = J dn k S(n • k)J • k8(n' ■ k)J' ■ k -> ^JJ'x 

J dX exp[iA(«; 2 - l)]w 2 dw [(n A n') 2 - 2iA/(L/c) 2 - A 2 /(L£;) 4 ] ~ 1/2 

In the limit L — > oo it results 

J = JJ7|nAn'| 

which can be obtained in a simpler way. 

Now we must integrate over r and r', times an oscillating factor 
e io; fc (r-T ) j n conditions that have been chosen the motion takes 
place in a plane, so n is characterized by a unique angle ip and n' by 
ip' hence the integration over time amounts at an angular integration, in 
fact (nAn') 2 = (sm(ip — ip')) 2 and moreover J = ip and J' = ip'. So 
we must integrate X in dip , dip' from to some final angle Defin- 
ing 7 = ip — ip 1 we see that the integrand shows, in the limit L — > oo, 
a singularity for 7 = 0, so we perform the integration from to |7r 

because the domain which do not include zero has no singular behavior, 
the oscillating factor is approximated with its value on the singular point 
r = r', so that the exponential factor reduces to 1. Then the integral is a 
complete elliptic integral, [12] which can be conveniently expressed in term 
of the hypergeometric function. In fact the integration in d^ is: 

x^JT^Tq Wlth ° = -(lSF-W (A3) 

The result of the integration is 



^"WoTT 2Fl (^ ;1; gTi) 

and in the limit L — > 00, which gives Q — > we get; 

X 7 -> ^[2 In 2 + In Life] . 

Since the dominant term in the limit is independent of the auxiliary pa- 
rameter A the rest of the integration in eq. (A.2) is straightforward and it 
gives 

l=VV2\n4kL (AA) 

A comment: what is excluded is the possibility that the magnetic field 
should perform more than a complete rotation, this would destroy the 
correspondence ip = ip' «-> r = r' . 
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